CALDERON COMMUTATORS AND THE CAUCHY INTEGRAL ON LIPSCHITZ 
CURVES REVISITED III. POLYDISC EXTENSIONS 
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Abstract. This article is the last in a series of three papers, whose scope is to give new proofs 
to the well known theorems of Calderon, Coifman, Mcintosh and Meyer [1], [6], [7]. Here we 
extend the results of the previous two papers to the polydisc setting. In particular, we solve 
completely an open question of Coifman from the early eighties. 



1. Introduction 

The present article is a natural continuation of the previous [19], [20] and is the last paper 
in the sequel. The goal of it is to show that the method developed in these papers to give 
new proofs to the LP boundedness of the Calderon commutators and the Cauchy integral on 
Lipschitz curves [1], [6], [7], can be used to extend these classical results to the n-parameter 
polydisc setting, for any n > 2. 

Suppose that F is an analytic function on a disc of a certain radius centered at the origin in 
the complex plane and A a complex valued function in R", so that a ^' A dx e L°°(R") with an L°° 
norm strictly smaller than the radius of convergence of F. Define the linear operator C„ ^ by 
the formula 



:= p.v. I 

JR" 



Qiyl/O) '■= P- V - f( X + t) F 



— o ... o —A(x) 



dt\ dt n 



for functions of n variables f(x) for which the principal value integral exists, where A® denotes 
the finite difference operator at scale s in the direction of e,, given by 



AfB(x) :=B(x + sed - B(x) 

and e\, e n is the standard basis in R". 

The main theorem we are going to prove is the following. 

Theorem 1.1. The operator C n j^ extends naturally as a bounded linear operator from L P (R") 
into L p (R n )for every 1 < p < oo. 

This answers completely an open question of Coifman from the early eigthties [10], [1 1]. The 
case when the L°° norm of „ d " \ is small and the generic n = 2 case, have been understood 

dx[...dx„ e> ' 

earlier by Journe in [10] and [11] respectively. Our proof is quite different from the approach 
in [10], [11] and works equally well in all dimensions. In fact, as we will describe in the last 
section of the paper, much more can be proved in the same way. Not only the operators of (1) 
are bounded, but also (for instance) those given by expressions of type 

l 
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f^p.v. f(x + t + s)F 



A? Ag> Ag> 

o o o A(x) 

h h s x s 2 



dt\ dt 2 ds\ ds2 
t\ h s\ s 2 



1^±1H11± (2) 



d*A 

dx^dx^ 



and their natural generalizations. Of course, in (2) one has to assume this time that 

L°°(R 2 ). When F(z) = z d with d > 1 the operator in (1) is the natural n-parameter extension of 
the Jth Calderon commutator, whereas for F(z) = one obtains the n-parameter generaliza- 
tion of the Cauchy integral on Lipschitz curves [1], [6], [7]. 

For simplicity, we shall denote from now on with C n ^ A the n-parameter Jth Calderon com- 
mutator. It is easy to observe that when fix) and A(x) are particularly given by 



and 



one has 



fix) = fl(Xi) ■ ... ■ f n (x n ) 



A(x) =A l (x 1 ) ■ ... -A n (x n ) 



C n ,d,Af( x ) ~ Cl,d,Aif( x " ••• ' ChdAnfn( x n)- 

To motivate the introduction of the operators C n> p,A one just has to recall the context in which 
the original Calderon commutators appeared [1], [2], [6]. If one tries to extend Calderon's 
algebra to R" and to include in it pseudodifferential operators containing partial derivatives, 
one is naturally led to the study of the operators in (1) and their natural generalizations. 

It is clear and very well known that to prove statements such as the one in Theorem 1.1, one 
needs to prove polynomial bounds for the corresponding Calderon commutators C n ^ A . More 
specifically, Theorem 1 . 1 reduces to the estimate 



\\C n ,d,Af\\p<C{n,d)-C{p) 



d n A 



(3) 



dx\... dx n 

for any / e LP, where C(n, d) grows at most polynomially in d l . 

The argument of [10] to prove the small L°° norm theorem used an induction on the dimen- 
sion n. We work instead directly in R" and since our method is essentially similar in every 
dimension, to keep the technicalities to a minimum, we chose for the reader's convenience to 
describe the proof of the main Theorem 1.1 in the particular case of the plane R 2 . However, it 
will be clear that the same proof works equally well in every dimension. 

So from now on n = 2 and the goal is to prove the corresponding (3). The operators C 2 ,d, A 
that we would like to understand, are given by 



'24 ,A. 



f(x) = p.v. J^ 2 f(x + t) 



'A (1) . 
=^ o ^Aix) 



V 



A? 
ti 



dt\ dt 2 
h h ' 



(4) 



lr This reduction is a simple consequence of the fact that if one writes the analytic function F as a power series, 
the generic operator C„,f,a itself becomes a series involving all the commutators C n ,d,A- The polynomial bounds 
are necessary for this series to be absolutely convergent. 
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a2 a 

If a := ° t , then one observes that 

v(D A (2) 



Aft A ?2 f 

h h J[o,i] 2 



(5) 



As in [20], using (5) d times, one can see that if a and / are Schwartz functions, the implicit 
limit in (4) exists and can be rewritten as 



Jr m+2 



where 



m 2M (^,^,...,^ d ,j],ri u ...,j] d ) := m hd (^,^ u ...,£j d ) ■ m w (i/,i7i, ...,%) (7) 
with m w (^, fx, &) and m M (?7, 771 , i^) given by 



I sgn( 

J[0,l] rf 



i(f + aifi + ... + a d £ d )dai...da d 
and 

I sgn(?7 +^i?7i + ... +/3 d j] d )d/3 1 ...d/3 d 

J[0,l] d 

respectively. Because of the formula (6) C 2 ,</ can be seen as a (J+ l)-linear operator. However, it 
is important to realize (as in [20]) that even though its symbol m 2 ^ has the nice product structure 
in (7), it is not a classical bi-parameter symbol, since m xd itself is not a classical Marcinkiewicz 
Hormander Mihlin multiplier 2 . As a consequence of this fact, the general polydisc Coifman 
Meyer theorem proved in [21], [22] cannot be applied in this case. The strategy would be to 
combine the techniques of [21], [22] with the new ideas of [19], [20] and to show that (together 
with some other logarithmical estimates that will be proved in this paper) they are enough to 
obtain the polynomial bounds of (3). Given these remarks, it would clearly be of great help for 
the reader, to be already familiar with our earlier arguments in [19], [20]. 
We will prove the following 

Theorem 1.2. Let 1 < p\,...,p d +\ < oo and 1 < p < oo be so that \/p\ + ... + l/p d+ \ = l/p. 
Denote by I the number of indices ifor which p t ± oo. The operator C 2 , d extends naturally as a 
(d + X)-linear operator bounded from L Pl x ... x L Pd+l — > LP with an operatorial bound of type 

C(d)-C(l)-C( Pl )-...-C( Pd+1 ) (8) 
where C(d) grows at most polynomially in d and C(p,) = 1 as long as p { = oofor 1 < / < d + 1. 



m\4 is of course the symbol of the one dimensional dth Calderon commutator [20]. 
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The above Theorem 1.2 is the bi-parameter extension of the corresponding Theorem 1.1 in 
[20]. If we assume it for a moment, we see that (3) follows from it by taking p x = p and 
Pi = ... = Pd+i = °°- 

To show Theorem 1.2 we will prove that for every 1 < i < d + 2 and for every 0i, (f> d+ \ 
Schwartz functions, one has 

||Ci(0!, <^ +1 )|| p < C(d) • C(Z) • C(p0 • ... • C(p d+1 ) • inu • ... • ||^ + ilU, (9) 

where {pj) d ^\ and /? are as before and (C^Sf are the adjoints of the multilinear operator C 2 ^. 3 
Standard density and duality arguments as in [20], allow then one to conclude that the estimates 
in (9) can be naturally extended to arbitrary products of U> and L°° spaces. 4 

Our plan for the rest of the paper is as follows. In the next section, Section 2, we describe 
some discrete model operators whose analysis will play an important role in understanding 
(9). In Section 3 we prove that the main estimates (9) can be reduced to a general theorem 
for the model operators. In Section 4 we prove the theorem for the discrete model operators of 
Section 2. In Section 5 we show logarithmical bounds for some shifted Hardy-Littlewood-Paley 
hybrid operators, which appear naturally in the study of the previous discrete models. Finally, 
in Section 6 we describe various generalizations of the main Theorem 1.1. 

Acknowledgements: The present work has been partially supported by the NSF. 

2. Discrete model operators 

As mentioned earlier, the main task here would be to describe some discrete model operators, 
whose analysis is deeply related to the analysis on (9). Because of the formula (6), we now know 
that 



C 24 = C U ®C U (10) 

and so one should not be at all surprised, to find out that these bi-parameter model operators that 
will be introduced, are in fact tensor products of the one-parameter discrete model operators of 
[20]. And also as in [20], these operators are not going to be (d + l)-linear, but /-linear instead, 
for some 1 < / < d + 1. The explanation for this is similar to the one in [20]. To be able to 
prove (9), one first decomposes C 2 ,d into polynomially (in d) many bi-parameter paraproduct 
like pieces and then estimate each such piece independently on d. To be able to achieve this, 
one has first to realize that one can estimate most of the L°° functions easily by their L°° norms 
and reduce (9) in this way to the corresponding estimate for some minimal /-linear operators. 
To prove the desired bounds for these minimal operators, one has to interpolate between some 
Banach and quasi-Banach estimates, as in [20]. The Banach estimates are easy, but the quasi- 
Banach estimates are hard. One has to discretize the operators carefully, in order to understand 
them completely. And this is (in a few words) how one arrives at the model operators. Their 
definition is as follows. 

A smooth function O(jc) of one variable is said to be a bump function adapted to a dyadic 
interval /, if and only if one has 

3 For symmetry, we aso use the notation C24 = C*^ 1 - 

4 The reader is also referred to our earlier [20] for an explanation of why does one need the larger set of estimates 
in (9) for C24 and its adjoints, even though one is interested in the more particular (3). 
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for all derivatives a satisfying \a\ < 5 and any large M > with the implict constants depending 
on it. Then, if 1 < q < oo, we say that \I\~ llq <S> is an L q normalized bump adapted to /. The 
function O(^) is said to be of *F type if Q>(x)dx = 0, otherwise is said to be of O type. 

A smooth function 0(x, y) of two variables is said to be a bump function adapted to the dyadic 
rectangle R = I x J if and only if it is of the form 0(x, y) = <&i(x) • <D 2 (y) with <&i(x) adapted to 
/ and <D 2 (y) adapted to J. If / is a dyadic interval and n an integer, we denote by I n := / + n\I\ 
the dyadic interval having the same length as / but sitting n units of length |/| away from it. 

Fix now 1 < / < d + 1 and ri! = (nj, rij), n/ = (nj, n^) arbitrary pairs of integers. Define 
also n /+1 := (0, 0). Consider families (O^ ) R for 1 < j < I + 1 of L 2 normalized bump functions 
adapted to dyadic rectangles R n = I n i x J n 2 where R = IxJ runs inside a given finite collection 

3 j i 

'R of dyadic rectangles in the plane. Assume also that at least two of the families (<E>^ ) 7 for 

J 

1 < j < I + 1 are of *P type and that the same is true for the families (${ )j for 1 < j < I + 1. 

; 

The discrete model operator associated to these families of functions is defined by 



Wi....,J» = Zij^^< 1 >"-^ < ^><- ( U) 



\R\V-W 
Ren 1 1 



The following theorem holds. 



Theorem 2.1. For any *mc/z a finite family of arbitrary dyadic rectangles, the l-linear operator 
T K maps L Pl x ... x L Pl — > L p boundedly, for any 1 < p ls < oo wzY/z l/^i + ...II pi = l/p 
and < p < oo, with a bound of type 



( 2 

O 



| 1 log 2 < n{ > • log 2 < > 



(12) 



where in general, < m > simply denotes 2 + \m\. And the implicit constants are allowed to 
depend on I. 

This theorem is the bi-parameter generalization of Theorem 3.1 in [20]. As pointed out 
there, standard arguments based on scale invariance and interpolation, allows one to reduce the 
above Theorem 2.1 to the more precise statement that for every fj e U> with \\fj\\ Pj = 1 and 
measurable set E c R 2 of measure 1, there exists a subset E' c E with \E'\ ~ \E\ so that 



Re<R 1 1 7=1 (13) 

where f M := xe- As in [20], the fact that one looses only logarithmical bounds in the above 
estimates, will be of a crucial importance later on. 



6 



CAMIL MUSCALU 



3. Reduction to the model operators 

The goal of this section is to show that indeed (9) can be reduced to Theorem 2. 1 or more pre- 
cisely to its weaker but more precise variant (13). In particular, one can find here a description 
of all the ideas that are necessary to understand why it is possible to estimate the biparameter 
Calderon commutators C 2 ,d with bounds that grow at most polynomially in d. 

The reader familiar with our previous work will realize that this section is in fact a tensor 
product of the corresponding section in [20] with itself. As there, the first task is to decompose 
C24 into polynomially many biparameter paraproduct like pieces which will be studied later 
on. 

Non-compact and compact Littlewood-Paley decompositions. Let Q>{x) be a Schwartz func- 
tion which is even, positive and satisfying <&{x)dx = 1. Define also ^(x) by 

= <D(x) - 1<D(|) 

and observe that y V(x)dx = 0. 

Then, as always, consider the functions ^(x) and O fc (x) defined by 2 ky V(2 k x) and 2 k <b(2 k x) 
respectively, for every integer k e Z. Notice also that all the L 1 norms of Q) k are equal to 1 . 
Since ^(x) = O^x) - ®k-i(x) one can see that 

k<k 

and so 

or equivalently 

keZj 

for almost every £ e R. On the other hans, as observed in [20], since ¥(0) 
write v F(^ r ) as 

= f 

for some other smooth and rapidly decaying function (p. 

These are what we called the non-compact (in frequency) Littlewood-Paley decompositions. 
The compact ones are obtained similarly, the only difference being that instead of considering 
the Schwartz function O before, one starts with another one having the property that suppO c 
[-1,1] and 0(0) = 1. 

As explained in [20], the advantage of the non-compact Littlewood-Paley projections is re- 
flected in the perfect estimate 



(14) 

= *F'(0) = one can 



|/*0*(jc)|<||/|| 



(15) 
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which plays an important role in the argument. 

The generic decomposition of C24. Using (6), if /, f\, f d+ i are all Schwartz functions, one 
can write the (d + 2)-linear form associated to C24 as 

I I sgn(£ + <x x %\ + ... + a d % d )da { ...da d \ • (16) 

irnjl+...+i; <i+ l=0 

I sgn(?7 +P1TH + ... + J3 d r ld )dj3 l ...dj3 d 

V[0,l] d 

/(£ if) f x (£ u m)-fd+\ (&+ WHd+i )d%dt;i . . .d% d+ 1 dr] dr\ x . . .dr] d+ 1 . 
Then, by using the Littlewood-Paley decompositions in (14) several times one can write 

1 = J] %({)%(ti)---%(td)% +l (t<i + 0. (17) 

io,/l,.. .,/<;+! eZ 

As in [20] since for every (d + 2) tuple (Z , l\, l d+ i) e Z d+2 one has that either / > l\, l d +\ 
or /1 > / , l d +i ••• or > Z , Zrf, fixing always the biggest parameter and summing over 
the rest of them, one can rewrite (17) as 

2/5 / (^)o / (f 1 )...o,(&)o / (^ + i)+ 
Z,0/(#O/(6)...O/(&)W& + i). 

Also as in [20], we use in (18) compact Littlewood-Paley decompositions for the £ and % d +\ 
variables and non-compact ones for the rest of them. Every single term in the decomposition 
(18) contains only one *P type of a function and we would like to have (at least) two. To be able 
to producte another one, one has to recall that £ + ^ + ... + g d+1 = 0. Taking this into account, 
let us take a look at the second (for instance) term in (18) in the particular case when 1 = 0. We 
rewrite is for simplicity as 

(19) 

We know from before that ¥(£1) = tffii^i) and so we can write 

= fi?(£i)(-£ - 6 - ... - = -fi^i) - 66%) - - - 

Using this in (19) allows one to decompose it as another sum of 0(d) terms, containing this 
time two function of ¥ type, since besides £i^(£i) one finds now either a factor of type £<E>(£) 
or of type f 7 7 (^ 7 ) for some j = 2, ...,d + 1. 

If one performs a similar decomposition for every scale Z e Z and each of the terms in (1 8) one 
obtains a splitting of the function l{f+f 1+ ...+f (/+1 =o} as a sum of 0(d 2 ) expressions whose generic 
inner terms contain two functions of *P type as desired. 

Since we are this time in the biparameter setting, one has to decompose l{,+, 1+ ... +%+1 =o} m 
a completely similar manner. Combining these two decompositions, allows us to rewrite the 
(d + 2) linear form of C24 as 
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X ( sgn(£ + audi + ... + a£ d )da x ...da d \ • (20) 

I sgn(77 +^i77! + ... + p d T] d )dPi... dfi d 

W[0,l] d 

<(^5?tfi)...5?(&)op F1 (fo,i)- 

^(^(^O.-.ojMop 1 ^!)- 

/(£ t/)/i (f i , 171 ). . ./d+i (&+ uld+i )d%d%i . . .d% d+ 1 ofy^ . . . d?7 d+ 1 , 
which completes our generic decomposition. 

Recall that at least two of the families (O^^y))^ for < j < d + 1 are of *P type and likewise 

at least two of the families (^J(i7j))k 2 for < j < d + 1 are of *F type as well. We denote those 
indices by i\,i 2 and 71,72 respectvely. There are several cases that one has to consider which 
correspond to the positions of these indices. We call an index intermediate if it is between 1 
and d and extremal if it is either or d + 1. In [19] we essentially witnessed two cases. Case 1 
was when at least one of the *F positions corresponded to an intermediate index and Case 2 was 
when both of the T 1 positions were extremal. Since we now work in the biparameter setting, 
there are as a consequence four possible cases of type Case i <S> Case j for 1 < i, j < 2. 

Case 1 <8> Case 1. Assume here that i\ = 71 = and i 2 = 72 = 1. As mentioned earlier, the fact 
that z'i = i 2 = is not important, they can be anywhere else in the interval [0, d + 1]. Also, the 
fact that the intermediate indices i 2 and 72 have been chosen to be equal is not important either, 
but we chose them so for the simplicity of the notation. As in [20] we would like now to expand 
the two implicit symbols in (20). 

As there, let us denote by £ := £ + a 2 ^ 2 + — + cc d ^ d and by rj := 77 + /3 2 i] 2 + ... + /3 d i] d and recall 
from [20] that the idea is to treat the first symbol of (20) as being dependent on the variables £1 
and £ and similarly the second symbol of (20) as being dependent on 771 and 77. Also, since most 
of our functions do not have compact support in frequency, we need to consider some other 
compact Littlewood-Paley decompositions. We first write as in [20] 

l= Z ? «)= Z -+Z- + z - 

h,h h«h h~h h»h 

which can be rewritten as 



Z O ri (^ n (^) + ^ 1) + Z ^nimr^i). (21) 

n n n 

Then, we consider an identical decomposition, but for the variables 77 and 77! this time, where 
the summation is indexed over the parameter r 2 . If we insert (21) into (20) it becomes a sum 
of three distinct expressions that generate the subcases l a , 1/, and l c respectively. If in addition 
one inserts the analogous formula of (21) for the variables rf and 77! into (20) as well, one ends 
up with nine biparameter subcases of type Case l a <g> Case l a , Case l a ® Case \ b and so on. 
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Case l a <8> Case l a . To analyze the impact that these extra decompositions have, we consider for 
simplicity (as in [20]) the particular term corresponding to k\ = k 2 = 0. However, the argument 
we use is scale invariant. 

Let us ignore the symbol in (20) for now and just concentrate on the remaining expression 
which becomes 



V n 



£ [o ri (£)¥ ri (£)] O^^O^ 1 ^)..^;^)^ 1 ^,!) ■ (22) 

Z--E- 

Vr 2 <0 r 2 >0 ) 

which allows us to split our existing subcase into four additional subcases. 

Case \' a <8> Case 1^. So this corresponds to the situation when both r\ and r 2 are negative. As in 

[20], using the fact that ^(£1) is compactly supported and given that <S>Q l (gi) is also of ¥ type 
(in fact it is of the form one can rewrite the i; part of (22) as 

£ 2 ri O ri (|)Op(^(^ 1 )...Op(^)6p ri (^ +1 ) = 

n<o 

J> [o ri (J)f^(^)] • 5p(^(^ 1 )...6p 1 (^ +1 )O ri (?) 



for naturally chosen compactly supported functions O ri (£), W r [ (£1) and 
This allows us to split the symbol 



I sgn( 
Jo 



as a double Fourier series of the form 



y a 

/ 1 n 



"i 



(23) 



where the Fourier coefficients satisfy the quadratic estimates 

1 1 



C- 1 - 



= O- < 



< n > z < > w 

for an arbitrarily large number # > 0. See [20] for these important estimates. 



(24) 
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Clearly, there are similar calculations that one can make for the r\ part of (22). Using both of 
them, one can see that the particular contribution of \' a <g> Y a in (20) (at scale 1) becomes 

f f y r i y y a. y cj = f 



[o ri (^)-o r2 (^] 



(25) 



/(£, t])M^u m)- ■fd+i(€d+und+i)dt;dt;i... d^ d+ idr]drii. ..di] d+1 da 2 ...da d d/3 2 ...dfi d . 
Now, if one fixes a,/3, r l5 r 2 ,7z, tTJ, the corresponding inner expresion in (25) becomes 



f 

;;+r;i+...+iy +1 =0 



(26) 



'71 



<D ri (£ + a 2 £ 2 + ... + ad g d ) ■ ® r2 (n +J3 2 T] 2 + ... +Md) 



dgd£i...d£ d+l dr]dT] l ...dT] d+ i. 
To be able to continue the calculations we need the following lemma. 
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Lemma 3.1. If F, F\, F d+i , <D, O are all Schwartz functions, then one has 




F{%, Tj)F 1 (^ l ,T] l )...F d+l ^ d+l ,T] d+1 ) 





(£>(a£ + a\i;\ + ... + a d+l ^ d+ {) ■ d?(br] + b x r\ x + ... + fi d+ ir] d+1 )d^ 1 ...d^ d+ idr]dT] 1 ...dr] d+1 = 



a F(x l -ati,X2-at 2 )F l (xi-aiti,X2-bit 2 )...F d+l (x 1 -a d+ iti, X2-b d+l t2)®(h)Q>(t2)dx l dx 2 dtidt2 

Jr 4 

for a, a\, a d +i,b, b\, b d+ \ arbitrary real numbers. 

This Lemma 3 . 1 is the biparameter extension of Lemma 4. 1 in [20] and since its proof doesn't 
use any new ideas it is left to the reader. As pointed out in [20] there is also a natural general- 
ization of it, which states that the formula works for more than two averages (so one can take 

an arbitrary number of O functions and another arbitrary number of <D ones). 

As in [20], if G is now an arbitrary Schwartz function and a a real number we denote by G a 
the function defined by 



Alternatively, one has G a (x) = G(x - a). Using the above lemma and this notation the previous 
(26) can be rewritten as 






(f d+l *® l ' d+l ®®l' d+l )(x l ,X2> 



Q? n (ti)<b r2 (t2)dt i dt2dx l dx2 = 
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®o(ti)Qo(t2)dtidt2dxidx2 



10^ 20 2=^\ 



(27) 



7=2 ^ ' 



<Do(/i )Oo(? 2 )<i?i dtqdx\ dx2 ■ 
We have to remember now that all the calculations so far have been made under the assumption 
that k\ = k 2 = 0, but they can clearly be performed in general and then the analogous formula 
of (27) is 



1,0,- 



2,0, 



n ~*2 



\f*%/ 2Kl+n ®% 2 '^ 



(Xi,X 2 ) 



(28) 



2,1, 



;*2 +r 2 



(xi,x 2 )- 



11 /y*<^ a '®^ W7 K-v..-v2)- 



*2 



7=2 



<^o(h)^>o(h)dtidt2dxidx2- 

In conclusion, if one denotes by a = (a 2 , a d ) and by = (/3 2 , ...,/3 rf ) one sees that the part of 
(20) that corresponds to Case \' a ® Case Y a can be written as 



| J 2 r < \^ 2 n C r[ ^ C n Q r m,n,ni,n,m4$,h,t2 



(29) 



where c 2 1 ^ 2 '"'" 1 '"'" 1 ' a!,p '' 1 '' 2 i s t h e operator whose (d + 2) linear form is given by summing over 



K n<0 r 2 <0 n,n\ 

<&o(h)<&Q(t 2 )dhdt2dad$ 

Vi ,r2,n,nl ,7j,n] ,c?,/J,fi ,?2 
'24 

k\,k 2 the inner expressions of (28). 

To prove (9) for (29) one would need to prove it for the operators c^'J 2 '"'" 1 '"'" 1 '"^' 11 ' 12 with 

bounds that are summable over r { , r 2 ,7T, n{,n, n x and integrable over a, ft, h,t 2 . It is clear that 
these operators are essentially biparameter paraproducts and therefore one expects that the 
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method of [21], [22] should be used. That indeed will be the case, but on the other hand the 
appearence of all these parameters mentioned earlier, have the role to shift the implicit bump 
functions which appear in their definitions and as a consequence this time one has to be very 
precise, when evaluates the size of their boundedness constants. 

Since in our case away from the indices and 1 all the bump functions are of <D type, the idea 
is to use the perfect estimate (15) (or rather, its biparameter variant) to bound all the functions 
which are in L°°. To be more specific, let us denote as in [20] by S the set of all the indices 
2 < j < d for which pj ± oo. Set / := |S| + 2 and freeze all the L°° normalized Schwartz 
functions f corresponding to the indices in {2, d} \ S . The new resulted operator is a minimal 

/-linear operator which will be denoted by Q l ^ r ^ n ^ n ^ a $^^ _ 

Shifted hybrid maximal and square functions. It is now time to recall a few basic facts 
about biparameter paraproducts, to be able to go further. Consider two generic families of L 1 
normalized bump functions (<& J kl ) kl and (Q> J k )k 2 for 1 < y < / + 1 so that in each of them for two 
indices the corresponding sequences are of *P type. 

An /-linear biparameter paraproduct is an /-linear operator whose (/ + l)-linear form is given 

by 

Let us first assume that we are in a case similar to the one considered before and that the ¥ 
functions appear for the indices j = 1,2. Then, one can estimate the absolute value of (30) by 



l+l 



SSifflxuXz) ■ SS(f 2 )(x 1 ,x 2 ) ■ Y\ MM(f j )(x 1 ,x 2 )dx 1 dx 2 
where in general MM(f)(xi,x 2 ) and SS(f)(xi, x 2 ) are defined by 



MM(f)(xi, x 2 ) = sup \f * O fcl <g> ® k2 (x u x 2 )\ 

k\ ,k 2 

and 



SS(f)(x u x 2 ) = 



(3D 



1/2 

Yj\f*®h®®k 2 (Xl,X 2 tf 

\klM2 

respectively. In order for (31) to make sense, we assume of course that both (O^)^ and (O fe ) fe 
there, are of ¥ type. Since both MM and SS are known to be bounded in every LP space for 
1 < p < 00, the above argument proves that our particular biparameter paraproduct in (30) is 
bounded from L Pl x ... x LP 1 into LP as long as l/pi + ... + l/pi = l/p and 1 < p\, pi, p < 00. 

As one can imagine, the above / 2 x/ 2 x/°° argument can be twisted, in which case one naturally 
obtains hybrid maximal and square functions of type MS and SM defined by 



MS(f)(x u x 2 ) = sup 

h 



( V /2 
*<!>*, ®^ k2 {x u x 2 )f 

V k 2 ) 



(32) 
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and 



SM(f)(x u x 2 ) = 



J]sup \f*O kl ®a> k2 (x u x2)\ 2 



1/2 



(33) 



respectively. One has to assume in (32) that the family (O fe ) fe is of *F type and that (O kl ) ki is of 
*F type in (33), for both expressions to make sense. 

As observed in [22] all these hybrid operators are bounded in LP for 1 < p < oo as well 
and as a consequence, one can bound every biparameter paraproduct in arbitrary products of LP 
spaces, as long as all of their indices are strictly between 1 and oo. 

This discussion shows that in order to understand the operator Q l ^ r2,n,nx,n ' nx,a ^ tutl (and of 
course, all the other possible ones) one has to understand how to bound not only the above 
operators, but also their shifted analogs of type M ni M ni , S ni S n2 , M ni S„ 2 and S ni M ni which are 

defined similarly, but with respect to the shifted functions (O^* 1 ) kl and (<D* \ 2 . 

In [19] we understood completely the one-parameter shifted maximal and square functions 
M n and S n and proved their boundedness on LP spaces with operatorial bounds of type 0(log < 

n>)\ 

Now the arguments of [19] and [22] show that their hybrid biparameter analogs mentioned 
before, will also be bounded on LP spaces with operatorial bounds of type 0(log 2 < ri\ > log 2 < 
n 2 >), as long as one can prove logarithmic bounds for the so called Fefferman-Stein inequality, 
namely 



( N 


1/2 




( N 


1/2 


X \M n fj\ 2 




< C p log 2 < n > • 






\j=l J 




p 


V;=l ) 





which should hold true for every 1 < p < oo. 

This inequality will be proven in detail in a later section. Until then, we will use freely all 
these logarithmic bounds. 

Banach estimates for Q l ^ r 2,n,n x ,n,n x ,a$,t u t 2 ^ Q| ven ^ logarithmic bounds for the shifted maxi- 
mal and square functions described earlier, it is not difficult to see (as in [20]) that the operator 
^un,n,n u n,ni,a$,ti,t2 ^ s m( j ee( j bounded from L Si x ... x L s ' into L s as long as 1 1 s\ + ... + 1 /si = 1/s 
and 1 < s\, si,s<oo with operatorial bounds no greater than 

/ =; =; \2/ 

(C < r x >< r 2 > log < n > log < n > log < n\ > log < m > log < [ti] > log < [t 2 ] >) . 

(35) 

And this contribution is perfect, given the extra factors 2 n , 2 n that appeared before (recall 
that both ri,r 2 are negative in our case) and the quadratic decay in n, n^n, h~y. 



'The logarithmical bounds for M n are due to Stein and can be found in [24] Chapter II. 
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Quasi-Banach estimates for Q l ^ ri,n,nx " n,nx,a ^ Ml . Assume now that the index s above satisfies 
< s < oo and so it can be sub-unitary. We would like to estimate the boundedness constants 
of 



C 



•I,r\,r2,n,n\,n,n\,&fi,t\,t2 . 



2,d 



: If 1 x ... x V 



(36) 



This time one has to discretize the operators in the X\,X2 variables and then take advantage 
of the general result in Theorem 2.1. Arguing as in [20] we see that the problem reduces to 
estimating expressions of type 



1 1 Y 1 

2^2^ £-i |/?|('-i)/2 



n 

PS 



</,-, G>; j ® o 2 ; 7 



\(f d+ u<S>]' d+1 ®<Z> 2 / +1 )\ 



where the sum runs over dyadic rectangles of the formi? = / x J. By applying Theorem 2.1 we 
see that the operatorial norms of (36) can be majorized by 



(l 6n 2 6ri log < n > log < n > log <n x > log <n x > log < [t{\ > log < [t 2 ] >) 

and the same is true for all its adjoint operators. In the end, by using the same interpolation 
argument as in [20], one can see that the operator Q l sy^ n > n ^ n > n ^ a $> t ^ satisfies the inequality (9) 
with bounds that are clearly acceptable in (29) as desired. 

These complete the discussion of Case Y a ® Case Y a . The rest of the cases can be treated 
similarly after certain adjustments. Since all of these adjustments have been described carefully 
in [20], the only thing that is left is to realize that they work equally well in our tensor product 
framework. The straightforward (but quite delicate) details are left to the reader. 



4. Proof of Theorem 2.1 

The proof of Theorem 2.1 is based on the method developed in [21] and [22]. First, we need 
to recall the following lemma whose detailed proof can be found in [22]. 

Lemma 4.1. Let J c R be an arbitrary interval. Then, every bump function <pj adapted to J 
can be written as 



= J] 2- 1000, ¥, (37) 

where for each k e N, 0* is also a bump adapted to J but with the additional property that 
supp(cf)j) c 2 k J 6 . Moreover, if we assume JT (p j(x)dx = then all the functions (f>j can be 
chosen so that <pj(x)dx = Ofor every k e N. 



'2 k J is the interval having the same center as J and whose length is 2 k \J\. 
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Fix now the normalized functions f, f and the set E as in (13). Using the above lemma, 
one can estimate the (/ + 1)- linear form on the left hand side of (13) as 

|a*(/i,...,/m)I < y 2- 1000 ^' V -l—\{f u o l R )\...\{f h <nL >IK/, + i,< a >l 

&N* ™ lRl (38) 

where the new functions O^ 1 '^ have basically the same structure as the old O^ 1 but they also 
have the additional property that supp($£ u ) c 2 k R. We denoted by 2 l R := 2 kl I x 2 k2 J, Jc = 

(h,k 2 ) and \k\ = k\+k 2 . 

As before, the form (38) will be majorized later on by tensorizing two separate I 2 x I 2 x 
/ M x ... x r estimates with respect to parameters / and J. As a consequence, for every index 
1 < j < I + 1 there are hybrid square and maximal functions naturally attached to that position 
which we denote by (M - S)j. More specifically (M - S)j can be either the discrete variant 
of M n iM n i or 5 n i5 n 2 or M n i5 n 2 or 5 n iM n 2 depending on the positions of the corresponding ¥ 
functions. For simplicity, we do not write explicitly the dependence of these functions (M -S ) 7 
on the shifting parameters n 7 . Recall also that each of them comes with a boundedness constant 
which is no greater than 0(log 2 < nj > log 2 < n 2 >) 7 . 

We construct now an exceptional set as follows. For each k eti define 



n -sfl = \J{(*>y) e R 2 : -S)j(fj)(x,y) > C2 5 ^ log 2 < nj > log 2 < n 2 >}. (39) 
Also, define 



= {(x,y) e R 2 : MM(x a _J(x,y) > ^} (40) 



and then 



Cl_ m = {(x,y) e R 2 : MM(xa_ 5 ^(x,y) > ^}. (41) 

Finally, we denote by 

n = U 

It is clear that \Q.\ < 1/2 if C is a big enough constant, which we fix from now on. Then, define 
E' := E \ Q. and observe that \E'\ ~ 1. 

Fix then kN and look at the corresponding inner sum in (38). We split it into two parts as 
follows. Part I sums over those rectangles R with the property that 

*nfi^*0 (42) 

while Part II sums over those rectangles with the property that 

7 It is a standard fact that the continuous and the discrete variants of these operators behave similarly. 
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R n = 0. (43) 

We observe that Part II is identically equal to zero, because if R n ^ then R c Q_ 5 ||| 

and in particular this implies that 2*7? c Q_ 5 ^ which is a set disjoint from It is therefore 
enough to estimate Part I only. 

Since R n Cl c , * 0, it follows that < ± or equivalently, |tf n C2 C J > ^r|i?|. 

-5|*| |R| 2/ 1 _ 5 | fc | 1 2/ 

We are now going to describe / + 1 decomposition procedures, one for each function fj for 
1 < j < I + 1 . Later on, we will combine them, in order to estimate our sum. 
Independently, for every 1 < j < I, define 

C2 5 ^ log 2 < n 1 > log 2 < n 2 > 
= ^ R 2 : (M - S)j(fj)(x, y) > L >—} 



-5|*;|+1 

and set 



then define 



K 3 . = {R e <R : \R n Q y . | > 

-5W+1 -5W+1 1 2/ 



C2 51 ^ log 2 < n 1 > log 2 < n 2 > 

fi -5i*i +2 = {(x > y) e r2 : (M -sWiXx-y) > ^ — l 

and set 

ft j „ = {Re ( R\'R j . :|/?nn y '. | > 

-5|*|+2 -5\k\+l ~5\k\+2 21 

and so on. The constant C > is the one in the definition of the set E' above. Since there are 
finitely many rectangles, this algorithm ends after a while, producing the sets {Q, J S .} and {^.} 
such that 7? = \J SJ K SJ . 

We would clearly like to have such a decomposition available for the last function fi+\ as 
well. To do this, we first need to construct the analogue of the set ^_ 5 ,a, for it. Pick N > a big 
enough integer such that for every R e 7? we have \R n Q^ c | > ^-\R\ where we defined 

a l l l N = {(x,y) e R 2 : (M - S) l+l (f l+l )(x,y) > C2 N }. 
Then, similarly to the previous algorithms, we define 



and set 



then define 



and set 



C2 N 

= {(x,y) e R 2 : (M - S) M (f M )(x,y) > — } 



<R l ^ +l = {Re<R:\Rna l _ + ^ l \>-\R\}, 



C2 N 

Q% 2 = (ieR 2 :(M - S) l+l (f M )(x,y)) > 



18 



CAMIL MUSCALU 



<R!: l N+2 = {ReK\ <R M N+l : \R n Q^ +2 | > ^\R\}, 

and so on, constructing the sets {O. 1 ^ } and } such that K = U^K'** . 
Then we write Part I as 

Sl,...,S;>-5|£|,S/+i>-iV Re ^si,».,s;+i 

where ^. S i,...,. V;+1 := ^ n ... n Now, if 7? belongs to 7? Sl ,...,. v , +1 this means in particular 

that R has not been selected at either of the previous sj - 1 steps for 1 < j < I + 1, which 
means that \R n flj^l < \R n Q'+'.J < or equivalently |/? n ft^l > 2bl|J?|, 

|# n > ^1^1 • But this implies that 

|7?nQ^_ 1 n...nQ^: 1 |>i|7?|. (45) 

In particular, using (45), the term in (44) is smaller than 

si,...,si>-5\lc\,si + i>-N * e ^ii ,...,*(+] 

Z L _ Z ^l</ 1 '< i >l"-l^<n ( >IK/-^< a >^^3')^ 

lT + ,n...n£2 , p c <r) i i 

Sl,...,S;>-5|*|,J/+i>-Ar s l+l- 1 Ke «si,..,s;+i 

X/+1 
fl(M-S)y(/ y )U,y)^ 

^....^-SW.^-iV i'-Kr^i W 
/ 

< Z l5lCkl 11 l0g2 < ^ > l0g2 < "J > r " " - ' 2 " V ' +I • I' 

si,...,s ;+ i>-5|£u + i>-W -/ =1 (46) 

where 



0^ := /?. 

On the other hand we can write 



< 1^1 < \{(*,y) e R 2 : MMCr n , )(x,y) > i}| 
< \a l si \ = \{(x,y) e R 2 : (M - SM/iXx, y) > ^|^}| < 2™. 



Similarly, we have 

for every 1 < j < / and also 
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S/+10! 



7+1 



for every a- > 1 . Here we used the fact that all the operators (M - S ) 7 are bounded on U as long 
as 1 < s < oo and also that — 1. In particular, it follows that 



|Ok f | < 2 ilPiei • ... • 2' we '2'" +iae ' +1 (47) 

for any < 9 X , 9 M < 1, such that Q x + ... + 6 M = 1. 
Now we split the sum in (46) into 



i 

£ 2 5 * 1 f] log 2 < n] > log 2 < n 2 > 2"« • ... • 2— | (48) 

il,...,J/>-5|l|,i; + i>0 -/' =1 

+ Z ^ II l0g2 < ^ > l0g2 < n ' > 2 " V1 " - " 2 " V ' +1 I D '*m i- 

il,...,J;>-5|l|,0>fj+i>-Ar ^' =1 

To estimate the terms in (48) we use the inequality (47) as follows. First, we choose 6\, ...,6i 
small enough so that 1 - pjOj > for every 1 < j < I. Because of this, 9 M can become quite 
close to 1. To estimate the first term in (48) we pick a very close to 1 so that 1 - a6i+\ > 0, 
while to estimate the second term we pick a large enough so that 1 - aQ M < 

With these choices, the sum in (48) is at most 0(2 mw n'=i log 2 < nj > log 2 < n 2 >) and 
this makes the expression in (38) to be 0(Fl'=i l°g 2 < n ] > log 2 < n 2 >) as desired, after 
summing over k e N . 

This ends our proof. 

5. Logarithmic al bounds for the shifted hybrid maximal and square functions 

To complete the proof of the main theorem, we need to demonstrate the logaritmical bounds 
that have been used for the shifted hybrid maximal and square functions. As we mentioned 
before, the arguments of [19] and [21] show that they would follow from the following loga- 
rithmical bound for the vector valued Fefferman-Stein inequality. 

Theorem 5.1. Let n e X be a fixed integer and denote by M„ the shifted maximal operator 
associated to n. Then, one has 

N N 

||(2 \M n fj\y /2 \\ p < C p log 2 < n > -||(2 \fj\ 2 ) ll \ (49) 

7=1 j=l 

for every N and any 1 < p < oo. 

Proof The proof is a combination of the classical argument of Fefferman and Stein [24] with 
the new ideas from [19]. A nice description of the Fefferman-Stein inequality is in Workman 
[26] and we follow that presentation closely. There are three cases. Clearly, n is supposed to 
be large, otherwise there is nothing to prove. Assume also that n is positive, since the negative 
case is completely similar. 
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Case 1: p = 2. This case is very simple and it follows immediately from the theorem in [19] 
which says that M n is bounded on L 2 (and in fact on any LP) with an operatorial bound of type 
0(log < n >). 

Case 2: p > 2. To understand this case one first needs to observe the following lemma. 
Lemma 5.2. The following inequality holds 



for every a > and measurable functions f and O, where M n has been defined to be 



k=0 

Proof To prove this lemma we need to recall a few facts from [19]. Denote by /„ maximal 
dyadic intervals selected with the property that 



Clearly, they are all disjoint and their union is equal to {x : Mf(x) > a}. Each /„ comes with 
[log 2 n] dyadic intervals of the same length attached to it, denoted by l\, /^ log2 " ] . More 
precisely, I k n lies 2 k steps of length \I n \ to the left of /„. It has been observed in [19] that 




(50) 



[log 2 n] 





(51) 



{x : M n f(x) >a}c[jl n ul [ n U ... U /) 



[log 2 n\ 



Using these, one can majorize the left hand side of (50) by 



a 






(52) 



Now, for every y e /„ one can see that 



Using this in (52) one immediately obtains the desired (50). 




The result of the above lemma implies that 
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M n : L l (R,M n ®dx) -> L 1,00 (R, \0\dx). 

Since O can be assumed not to be identically equal to zero, we know that Al n O > 0. In 
particular, we also have the trivial bound 

M n : L°°(R,M n ®dx) -> L°°(R, |<D|</jc) 
and by interpolation, we obtain the following L 2 estimate 

f \M n f(x)\ 2 Wx)\dx < f \f{x)\ 2 M n <b{x)dx. (53) 
JR JR 

Coming back to the proof of Case 2, since p > 2 we know that g := p/2 > 1. By picking an 

appropriate \\<!>\\ q > = 1 and relying on the previous (53), one can write 

n n r N 

||(£ \M n fj\ 2 ) l % = II J] \M n fj\ 2 \\ q = I |M„//)|0|^ 

7=1 7=1 JR 7=1 

„ iV N N 

^ I l//)M,<M* * II E ^l 2 ^ ' H Ml0|l «' ~ ll( Z l//) 1/2 l£ • IIMJI^w 

JR 7=1 7=1 7=1 

On the other hand, from the definition of M n and the result of [19], we know that 

[log 2 n] [log 2 n] 

UAU,w < J] l|M- 2 *ll ? w < Y log 2* < log 2 < n > 

/t=i /t=i 

which completes Case 2. 

Case 3: 1 < p < 2. The idea here is to prove the following end point case 



N N 



11(2 |M„//) 1/2 || 1)00 < log 2 < n > -||(2 l//) 1/2 ||i (54) 

7=1 7=1 

directly and then to apply standard vector valued interpolation with the corresponding L 2 esti- 
mate. 

To prove (54), let a > and denote by F(x) := (ZyLi l/y(^)l 2 ) 1/2 - Select maximal dyadic 
intervals /„ with the property 



-i- f F(x)dx 

Vn\ Jl„ 



> a. 



As before, we like to think of each /„ as being related to the dyadic interval /, having the same 
length as I n and lying n steps of length |/ n | to the left of it. If we denote by Q, := \Jj I n one has 
as usual 

ini = e i^i < - e r p ^ dx ^ -ii^ii- ( 55 ) 
j a j Ji„ a 

Observe that F < a on Q c and also that 
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a< hl„ 



F(x)dx < 2a 



because of the maximality of /„. 

Split now each f k as f k = f k + ft where f' k := fkXo? and ft := fkXo.. 

Contribution of [ft. One can write 

\{x : (J] WJ'jixtff 12 > a/2}\ < ^\\(^ \M n f;(x)\ 2 ) m \\ 2 2 

j j 

< — 2 log 2 < n > ||( V |^«| 2 ) 1/2 || 2 < log 2 <n>-A F 2 (x)dx < log 2 < n > -\\F\U 
a *-r* a Jqc a 

as desired. 

Contribution of {ft}. To estimate the corresponding contribution for {ft}, we have to be a bit 
more careful. Define first the functions gk by 

^ := E ( rn f (56) 

'-f \ l n\ Jl„ 



and after that G{x) : = (2j \gj(x)\ 2 ) 1/2 . Fix x e I n and observe that by Minkowski's inequality 
one can write 



G(x) = 



< 2a-. 



Using that G is supported in Q. and arguing as before, we have 

\{x : \M ngj (x)\ 2 ) 1/2 > a/2}\ < \ log 2 < n > I^WI 2 ) 1/2 II 



2 

j j 



= \log 2 <n> \\G\\\ < log 2 < n > |Q| < ^log 2 < n > \\F\\ X . 

a 1 a 

We would like now to compare M n ft(x) with M n g k {x) if possible. Denote by Q, the set 

h:={j3I n u3I l n U...U3lH og2 " ] 
i 

and observe that |£2| < log < n > |Q|. We will prove that for every x e Q c one has 

M n ft(x) < M ngk (x) (57) 

and this will clearly allow us to reduce the contribution of {ft} to the contribution of {gk} which 
has been understood earlier. Fix x e OF and x e J a dyadic interval, so that the corresponding 
rn f, \fl'(y)dy is different from zero. In particular, J n has to intersect Q. which is the support 
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of /,". Suppose now that / is so that J n n /„ 0. Then, one must have /„ c J n as the other 
alternative J n c /„ is not possible since x e Q. c . But his implies that 

-]- f = -i- r i g ,(x)i^ 

which is enough to guarantee (57) and end our proof. 



6. Generalizations 

The goal of this section is to point out that virtually all the earlier generalizations that we 
described in the previous [19] and [20], have natural extensions in this multi-parameter world 
and can be proved by the same method. We give here just two samples and leave the rest (and 
the straightforward details) to the imaginative reader. Suppose for simplicity that we are in R 2 
and denote by D\ := ^ and D 2 := A direct computation shows that the double commutator 
[\D 2 \, [|Z)i |, A]] can be rewritten as 



[\D 2 \,[\D 1 \,AWix) = p.v. J(x + t) -^o-^A(x) 

JR 2 ^ *i h 



dty dt 2 



(58) 



which is precisely the bidisc extension of the first commutator of Calderon. There is of course 
a similar formula available in every dimension. 

Theorem 6.1. Let a\, a n be real numbers, all different from zero. The expression 



p.v 



fix + 1) 

d"A 



( a (1) a (n) 
O ... O A(x) 



\ 



dt\ dt n 



viewed as a bilinear map in f and d ^' A dx is bounded from L p xL q into U for every 1 < p, q < oo 
with l/p + l/q = \/r and 1/2 < r < oo. 



The particular case q = oo is in Journe [10] but the rest of the estimates seem to be new. 
Then, one can also observe by a direct calculation that 



=p - v -L fix 



[|D 2 |,[|D 1 |,[|D 2 |,[|D 1 |,A]]]]/(oc) 

A (2) a(D a (2) 



+ 1 + s) 



{ ^ A- A- A, 
o o o A(x) 



dt\ dt 2 ds\ ds 2 
h h s\ s 2 



{ h h s v s 2 

which is the bidisc analogue of an operator introduced in [20]. As we promised, we record now 
the following theorem. 

Theorem 6.2. Let F be an analytic function on a disc of a certain radius centered at the origin 
in the complex plane and A a complex valued function in R 2 so that 6 L°°(R 2 ) with an L°° 
norm strictly smaller than the radius of convergence of F. Then, the linear operator 

((2) a(1) a (2) \ 



f 



P.v. ^ fix 



+ t+s)F 



o o — -!- o -^A(x) 
h h si s 2 



dt\ dt 2 ds\ ds 2 
h t 2 si s 2 
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is bounded on LP (R 2 ) for every 1 < p < oo. 
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